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ABSTRACT 

We reconsider the origin and processing of dust in elliptical galaxies. We theoretically formu¬ 
late the evolution of grain size distribution, taking into account dust supply from asymptotic 
giant branch (AGB) stars and dust destruction by sputtering in the hot interstellar medium 
(ISM), whose temperature evolution is treated by including two cooling paths: gas emission 
and dust emission (i.e. gas cooling and dust cooling). With our new full treatment of grain 
size distribution, we confirm that dust destruction by sputtering is too efficient to explain the 
observed dust abundance even if AGB stars continue to supply dust grains, and that, except 
for the case where the initial dust-to-gas ratio in the hot gas is as high as ~ 0.01, dust cooling 
is negligible compared with gas cooling. However, we show that, contrary to previous expec¬ 
tations, cooling does not help to protect the dust; rather, the sputtering efficiency is raised by 
the gas compression as a result of cooling. We additionally consider grain growth after the gas 
cools down. Dust growth by the accretion of gas-phase metals in the cold medium increase 
the dust-to-gas ratio up to ~ 10“^ if this process lasts > lO/iriu/lO^ cm“^) Myr, where tin is 
the number density of hydrogen nuclei. We show that the accretion of gas-phase metals is a 
viable mechanism of increasing the dust abundance in elliptical galaxies to a level consistent 
with observations, and that the steepness of observed extinction curves is better explained 
with grain growth by accretion. 

Key words: dust, extinction — galaxies: elliptical and lenticular, cD — galaxies: evolution 
— galaxies: ISM — methods: analytical 


1 INTRODUCTION 


In the nearby Universe, elliptical galaxies are known to have little 
ongoing star formation activity. However, elliptical galaxies are not 
completely devoid of g as: they are also known to have hot X-ray 
emitting halo gas (e.g. I^ Sullivan et al.ll200lh and cold gas (e.g. 
IWiklind. Combes. & Henk^ 1995h . Moreover, dust is detected for 
a large (>50 per cent) fraction of elliptical galaxies by optical ex¬ 
tinction (e.g. iGou dfrooii et alJll99 4al : lvan Dokkum & Franx|[T995l : 


[Ferrari et alJll99^: Tran et al. [2001 ) or far-infrared (FIR) emission 


(e.g. lKnapp et al.lll989t) . Dust mass is estimated to be lO'^-lO ^ Mq 
from the reddening in the optical jGoudfrooii et al.lll994ah . The 
dust mass estimated from FIR dust emission tends to be ev en larger 
(~ 10^-10^ M©; e.g. lUeeuw et al.ll200l : lTemi et al.ll2004ll . The ex¬ 
istence of dust is also important in determining the chemical prop¬ 
erties of the cold gas through dust-surface reactions and freeze- 
out of molecules onto the dust (Fabian, Johnstone, & Dainesll 19941 : 
IVoit & Donahueiri995h . In spite of such an important role of dust, 
the origin of dust in elliptical galaxies is still being debated, and 
there is no clear theoretical explanation for the total amount of dust 
there. 


* E-mail: hirashita@asiaa.sinica.edu.tw 


Because the stellar population is dominated by old stars whose 
ages are comparable to the cosmic age, the dust is predominantly 
supplied by asymptotic giant branch (AGB) stars in elliptical galax¬ 
ies. However, the supplied dust is quickly destroyed by sputtering 
in the X-ray-emitting hot gas. This destruction is so efficient that 
the observed dust mass cannot be explained by the dust abundance 
achieved by the bal ance between the supply from AGB stars and 
the destruction (e.g. jPatil et aDl2007h . Therefore, the existence of 
such an ‘excessive’ amount of dust has long been a mystery. 

Some authors argue that the dust existing in elliptical galax¬ 
ies is possibly injecte d from outside via the merging or accretion 
of ex ternal galaxies (Forbes! Il99ll : iTemi et aP |2004 iFuiita et ^ 
l2013h . The lack of correlation between dust FIR luminosity and 
stella r luminosity is also taken as evide nce of external origin of 
dust (Temi. Brighenti. & Mathewsll2007h : however, this argument 
implicitly assumes that the existing dust is tightly related to the 
stellar dust production activity. If dust is processed by mecha¬ 
nisms not related to stars, it may be natural that we do not find a 
correlation between dust emission and stellar luminosity. In par¬ 
ticular, dust may grow through the accretion of gas-phase met¬ 
als in the cold gas, and this dust growth is suggested to be the 
most eff icient mechanism of dust mass increase in nearby galax- 
ies (e.g. iPwe^i 19981 : lHirashitairi999l : IZhukovska. Gail. & Trieloffl 
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20081; Inoue 2011; Matts son & Andersenl 2012 

; Asano et alJ 

2013al; de Bennassuti et al.l2014l; Remv-Ruver et al. 

2014). Indeed. 

Fabian et al.l (Il994ll and IVoit & Donahue! (Il995l) 

pointed out a 


possibility of dust reformation in a c o oling flow. More recently, 
[Martini, Dicken, & Storchi-Bergmamil ( 12013h also suggested the 
importance of accretion in maintaining the dust abundance to a 
level consistent with observations by analogy with the Milky Way. 
Since the efhciency of dust growth is not necessarily related to the 
stellar properties of elliptical galaxies, the lack of correlation be¬ 
tween dust and stellar emissions does not necessarily mean the ex¬ 
ternal origin of dust. 

The interstellar dust is not only passively processed in the 
ISM, but could also actively regulate the physical state of the 
ISM, since dust is able to radiate away the energy obtained by 
collisions with gas particles. This is one of the most important 
cooling paths of the hot IS M, and referred to as dust cooling. 
[Mathews & Brighentil ( l2003h argued that dust cooling could oc¬ 
cur faster than the dynamical time in the galactic potential and 
thus possibly enable the dust to survive against sputtering. How¬ 
ever, they did not treat the dust supply from AGB stars consis¬ 
tently, and they assumed a relatively high initial dust-to-gas ra¬ 
tio (~0.01, which is comparable to the dust rich ISM seen in the 
Milky Way) in the hot gas. Since there has been some advance 
in the theoretical understanding of dust pro duction in AGB stars 
dperrarotti & Gaii2006l : IVentura et al.ll20l3) , it is worth reconsid¬ 
ering their model s using the new know l edge a bout the dust yield. 
Moreover, since [Mathews & Brighentil ( l2003h assumed a single 
grain size for dust cooling, it would be interesting to calculate the 
evolution of grain size distribution: indeed, the grain size distri¬ 
bution is of fundamental importance in comparison with observed 
extinction curves as we argue below. The rates of s puttering and 
dust cooling also depend on the grain size distribution dOwe kll987l: 
iNozawa. Kozasa, & Habel[2006l) . 

The evolution of dust is also link ed to active galact ic nucleus 
(AGN) activities in elliptical galaxies. |Temi et ^ 2007!) proposed 
the transportation of dust from dust reservoirs in the central re¬ 
gions to explain the extended dust emission in elliptical galaxies 
(see also [Mathews et ai]|2013h . In this context, the origin, produc¬ 
tion, and survival of dust are important in clarifying the origin of 
the dust reservoirs. As mentioned above, dust may also contribute 
to cooling. The cool gas that falls into the centre may fuel the 
AGN (e.g. IWemer et al.ll2014l) . which in turn works a s a heating 
source and reforms th e hot g as dCiotti & OstrikeilboOlh . Recently, 
IValentini & Brighentil d2015l) pointed out that the AGN feedback 
could also enhance cooling locally through compression. 

In this paper, we focus on the processing of dust in ellipti¬ 
cal galaxies, both in the hot and cold gas components. In order to 
avoid complexity arising from the AGN feedback, we do not treat 
the supply of hot gas through the feedback but start the calculation 
given the pre-existing hot gas. Instead, we do treat cooling of the 
gas in a consistent manner with the dust supply and destruction in 
the hot gas. Since dust cooling and dust destruction depend on the 
grain size distribution, we fully take the evolution of grain size dis¬ 
tribution into account. We additionally consider the impact of dust 
growth in the cooled gas not only on the dust amount, but also on 
the grain size distribution. Through the modeling of the above pro¬ 
cesses, we will be able to clarify the evolution and survival of dust 
in a cooling cycle of the ISM in elliptical galaxies. 

We emphasize that the treatment of grain size distribution is 
one of the most important features in our modeling. Indeed, there is 
an observational clue to the grain size distribution, which enables 
us to test our results. The wavelength dependence of extinction. 


the so-called extinction curve, is investigated for a statisti cal num¬ 
ber (>10) of elliptical galaxies: iGoudfrooii et al.l (Il994bll showed 
that the ratio of total to selective extinction Ry, which is a mea¬ 
sure of the flatness of optical extinction curve, tends to be smaller 
{Ry - 2.1-3.3) th an the mean value of the Milky Way extinc- 
tion curves (3.1; e.g. |Peil992l:IWeingart ner & Drainj200l l:lDrainel 


I 2 OO 3 I : [Fitzpatrick & Massall2007 ; Nozawa & Fukugita 2013ll . This 

implies that the typical grain size is small er in the elliptical galaxies 
than in the Milky Wav [Patil et al.l ( l2007h obtained a similar range 
{Ry - 2.03-3.46) for 11 elliptical galaxies. These extinction curves 
can be used to examine the evolution of grain size distributions 
through various processes in the ISM. 

The paper is organized as follows: we explain our models for 
the basic processes of dust evolution in the hot ISM in Section 
We show the results in Section [3] We additionally consider dust 
processing in the cold ISM, focusing on dust growth, in Section|4] 
In Section 0 some predicted features of the models such as dust 
abundance and extinction curves are compared with observations. 
Finally we conclude in Sectionj^ 


2 PROCESSING IN THE HOT ISM 

We explain the method of calculating the evolution of dust grain 
size distribution in the hot ISM in an elliptical galaxy. Our mod¬ 
els are composed of dust supply by AGB stars and dust destruction 
by sputtering. Cooling of the hot gas by dust thermal emission and 
gas line emission is also calculated in a consistent manner with the 
evolution of dust grain size distribution. To avoid the complexity 
arising from dynamical evolution determined by the gravitational 
field and the AGN heating, we only consider the evolution of a 
local ‘fluid element’ of the ISM without modeling the global dy¬ 
namics. The simplicity of this approach enables us to focus on dust 
processing which is determined by the local environment where the 
dust resides. 


2.1 Evolution of grain size distribution 


The grain size distribution at time t is defined so that n(a, f) da is 
the number density of grains in the range of grain radii between a 
and a-I-da. The time t is measured from the onset of dust processing 
in the hot gas (t is not the age of the galaxy). As mentioned above, 
we focus on a local region without considering the global gas and 
dust distributions within the galaxy. We also assume that the gas 
and dust are well coupled (i.e. once dust is injected into a certain 
gas element, that dust stays in the same element without any dis¬ 
placement; in other words, the Lagrangian trajectories of gas and 
dust are the same). Under these assumptions, the evolution of grain 
size distribution in the hot gas is described as 


9n(a, t) 
dt 


9 d In p„as 

H- [an(a, 1)] = S (a, t) + n(a, t) - 

oa dt 


( 1 ) 


where a = da/d? is the changing rate of grain radius by sputtering 
(d is negative; see Section l23t . 5 (a, ?) is the source term due to 
the dust production by AGB stars (Section|2j2]l, and Pgas is the gas 
mass density. The last term in equation Q is the increase of the 
grain number density by the change of the background gas density 
(AnnendixlAt: as formulated in Section lT4l we consider the change 
of gas density by cooling. We give the number density of hydrogen 
nuclei, n^, as an input parameter, and the relation between and 
Pgas is given by 


Pgas = pmunn. 


( 2 ) 
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Dust processing in elliptical galaxies 3 


where oth is the mass of hydrogen atom, and ji is the mean weight 
of the gas per hydrogen (we adopt yu = 1.4). 

Since we are often interested in the grain mass, we introduce 
the grain size distribution per unit grain mass hy defining h(m, t) Am 
as the number density of grains with mass between m and m + dm; 
i.e. n{m, t)Am = n(a, t]da. The grain mass is estimated as m - 
^na^s (i.e. grains are assumed to be spherical), where constant.? is 
the material density of the grain. For the calculation of grain size 
distribution, we adopt i = 3.3 g cm“^ based on a silicate material. 
Using th e fact that a is independent of a (Section [2.3b . we obtain 
(see also lHirashitall20 12 I) 


d(mn) 

At 


A(mh) 

Am 


1.-, p, , -dlnpga. 

—mn + mS (m, t) + mn -, 

3 df 


(3) 


where m = dm/At = 4na^sd and S (m, t) dm = S (a, t) da. 

The dust-to-gas ratio as a function of time, T>{t), is calculated 
by 


0(0 


1 

ftmiitiH 

1 

ftmuriH 



—na^sn{a, t)da 


mh(m, t) dm. 


(4) 


2.2 Supply from stars 

2.2.1 Formulation 


We consider the stellar population contributing to the dust supply in 
the local gas element considered above. Since we are interested in a 
system as old as the cosmic age, the current stellar population con¬ 
tributing to the dust formation is AGB stars, which originate from 
stars with m < 8 (m is t he mass at the zero-age main sequence) 
(e.g. lZhukovska et alfcOOSll . The dust mass density supplied from 
AGB stars per unit time, Pdust,AGB. is estimated as 


r^m 

Jmt. 


pdust,AGB = I md(m)(/i(m)p,(4ge - T,„)dm, (5) 

where is the turn-off mass (i.e. the mass of the star whose 
lifetime is equal to the galaxy age t^gP, m^p is the upper cut-off of 
stellar mass, m^(m) is the total dust mass produced by an AGB star 
with mass m, c/>(m) is the stellar initial mass function (IMF), p,(?age) 
is the star formation rate density as a function of galaxy age, and t„, 
is the lifetime of the star with mass m. We take the stellar lifetime 
from iRaiteri. Villata. & Navan3 ( Il996h . We adopt a Salpeter IMF 
[0(m) oc with m^p = 100 Mq and miow = 0.1 Mq, where 

OTiow is the lower cut-off of stellar mass, and the normalization of the 
IMF is determined by f m(p{m)dm = 1. The data of dust yield 
md(m) is described in Section [2.2.21 We assume that dust grains 
ejected from AGB stars are instantaneous ly (on a time-scale shorter 
than 10^ yr; Math ews & Brighenti|[i999l) injected into the hot gas 
dTemi et alj2007l) . We approximate the star formation history with 
an instantaneous burst at fage = 0, which is a good approximation 
for nearby elliptical galaxies: 


P*(tage) “ P*,tott^(tage). 


( 6 ) 


where p, tot is the total stellar mass density ever formed and S(t) is 
Dirac’s 6 function. Noting that only stars with mass m,^g^ are con¬ 
tributing to the current dust formation, we obtain from equations 
([3 and ([fill 


Pdust.AGB “ l^d(l^rage)^i^rage)P*.lot 


dm, 


At 


(7) 


Taking the mass loss of stars that have ended their lives into ac¬ 
count, the current stellar mass, p,, is estimated as 

p. = (l-'R)p,,o„ (8) 

where K is the returned fraction. We assume that K - 0.25 
dKuo et alj|20l^ (this value is derived f or ^ 1 Mpi, but i s 
not sensitive to fage as long as fage > 10^ vn lFlirashita & Kucl201 ih . 
Correspondingly, we adopt fage = 9.5 Gyr, the lifetime of m = 1 Mq 
star. We fix fage, since we are interested in the time evolution at 
f £ 10* yr, which is much shorter than fage. 

Now we formulate the source term in equation dU. This term 
expresses the dust input from AGB stars, and we hereafter refer to 
the dust supplied by AGB stars as AGB dust. We specify the grain 
size distribution of AGB dust in such a way that the source term in 
equation (ID is written as 


S (tZ, ?) — Pdustfd\ist{.^\ 


(9) 


where /dust (a) is the grain size distribution function whose normal¬ 
ization is determined by 


r 


-Ttasfia) da = 1. 


( 10 ) 


Or equivalently, it is possible to use equation © with S (m, t) = 
Pdust/dust(fff). where the normalization of mass distribution function 
fiusiim) is determined by mpjstCtn) dm = 1. For the grain size 
distribution, we consider two ca ses: one is a power-law form with 
the same power index as used in iTsai & Mathewa dl995ll . 


/dust(fl) = CmQ (0.001 qm < a < 0.25 qm), 
and the other is a lognormal form. 


_ ^ C, f [ln(a/ao)]2 

/dus.(«) = -exp- 


( 11 ) 


( 12 ) 


where Cm and Ci are the normalizing constants determined by 
equation Gil, and ao and cr are the central grain radius and the 
standard deviation of the lognormal distribution, respectively. The 
first one (e quation 1 11b is based on the one de rived for the Milky 
Way ISM jMathis. Rumpl. & Nordsiecklll977h . and is referred to 
as the MRN grain size distribution. The second one in equation 
Gi is based on the argument th at AGB stars produce large (0.1- 
1 qm) grains dAsano et alj2013 fl and reference s therein). We adopt 
the same parameters as in Asano et^ ( l2013bh : ao = 0.1 qm and 
cr = 0.47, which are based on t heoretical calculations of dust con¬ 
densation in AGB star winds bv lYasuda & Kozas3 ll2012[) . 


2.2.2 Dust yield of AGB stars 

We take mjfm) (the ejected dust mass by an AGB star for different 
progen itor mass) fro m theoretical calculation s bv izhukovska et aP 
(|20()8[| . who adopted [perrarotti & Gail! ( l2006h ’s framework to cal¬ 
culate the dust condensation in AGB star winds. Since we are inter¬ 
ested in nearby elliptical galaxies, it is reasonable to consider the 
solar metallicity case (Z = 0.02). The stellar mass around 1 Mq, 
which is contributing to the current dust production, is relevant to 
this paper. Their calculations indicate that around 4x lO^'^-lO^* Mq 
of dust is produced by a m = 1-2 Mq star. If we adopt Z = 0.015 
(0.03) instead, the dust mass becomes 3 times as small (large) as 
the case above. Therefore, we should keep a factor of 3 uncertainty 
in mind. I Ventura et al.l ( 12014|) also calculated theoretical dust yield 
for AGB stars, focusing on lower metallicities. Their highest metal¬ 
licity case (Z = 8 X 10“*) shows that the dust mass produced by an 
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4 Hirashita et al. 


AGB star with m = 1-2 Mq is 2 x 10 "*-10 ^ Mq, similar to that 
adopted above. 


2.3 Sputtering 


The changing rate of grain radius by sput tering, a, in equation ([T} is 
estimated by a fitting formula derived bv IXsai & MathewsI ( 1199^ : 


a = -h 



1 + 



(13) 


where the coefficients of this fitting function, h = 3.2 x 10“** cm"* 
s“* and T = 2 x 10® K, are determined to giv e a good fit to the 
sputte ring rate of both silicate and graphite in IPraine & Salpeten 
( I1979I) . This formula correctly reproduces the temperature depen¬ 
dence of sputtering in that it is efficient at T > 2 x 10® K and 
steeply declines at T < 10® K dNozawa et alj200^ . We also define 
the sputtering time-scale, Tsput = |a/d|, which is estimated aiT > T 
using equation 13 as 


”^sput 



yr. 


(14) 


This is consisten t with the values der ived by a detailed treatment of 
sputtering yield iNozawa et al .I2OO6I see their fig. 2). 


2.4 Cooling of the hot ISM 

Following iMathews & Brighentil l2003h . we assume that the gas 
pressure in the local region of interest is constant. The time evolu¬ 
tion ofthe_gas_tem£eratme7^rough isobaric cooling is described 
by iMathews & Brighentill2003h 


dr 2AM(r) 

df 5nk 


(15) 


where A,ot(r) is the cooling rate (energy lost per volume per time) 
including both gas and dust cooling, n is the number density of gas 
particles, and k is the Boltzmann constant. Assuming that H and He 
are fully ionized, we estimate that n = 2.3«h- 

The cooling r ate b y gas emission is taken from 
[Sutherland & Dotfitt] ll993h . Specifically, we adopt n,neAN(r) 
[AnIT) is the cooling function for gas cooling] for the cooling 
rate per volume with electron number density = 1.2nH and 
total number density of ions n, = 1.1«h- We stop the calculation 
when T reaches 10® K (this time is denoted as fs), below which 
the temperature evolution occurs much faster and the effect of 
recombination would change significantly. Since sputtering does 
not work at such a low temperature, the resulting grain abundance 
and size distribution are insensitive to the detailed treatment of 
thermal evolution at < 10® K. 

The gas also cools through ga s -grain collisions (e.g. 
Dalgarno & McCravl 1 19721: iDwelJ Il987l: ISeok. Koo. & Hirashital 
2015t) . The rate of this dust cooling is estimated as n^^nAdfr), 
where Ai(T) is the cooling fu nction for dust cooling. We compute 
Aj(r) based on iDwekI l ll987h under the dust abundance and grain 
size distribution calculated in this paper. We again adopt = 1.2nH 
and assume silicate as dust species (adop ting graphite instead does 
not significantly change the cooling rate: lDwe 

After all, we evaluate the total cooling rate as A,ot = -1- 

«(,«HAd. We assume a constant pressure, so that the density is si¬ 
multaneously varied as 


dlnpgas _ dlnT 
df df 


Table 1. Model parameters adopted from iTsai & Mathe^ il995li . 


Model 

(10*° Lq) 

P. 

(g cm“®) 

Pgas“ 

(gcm^®) 

T 

(10® K) 

a 

10.6 

1.19 X 10“®' 

1.43 X 10“®® 

0.922 

b 

3.31 

1.36 X 10“®° 

3.01 X 10“®® 

0.468 

c 

0.976 

1.83 X 10“*° 

8.74 X 10“®® 

0.252 


Note. The densities at the core radius are adopted. 

“The corresponding hydrogen number densities are nn = 0.0610, 0.128, and 
0.373 cm“® for Models a, b, and c, respectively. 

2.5 Models and initial conditions 

For simplicity, we only focus on a ‘typical’ single region in an el- 
lipti cal galaxy by adopting re presentative gas density and tempera¬ 
ture. iTsai & MathewsI lll99.5h considered three models of elliptical 
galaxies which are classified by 6-band luminosities, Lg. The mod¬ 
els are referred to as Models a, b, and c depending on the 6-band 
luminosity. Although we do not directly use 6 band luminosity, the 
quantities used in the models (p.,Pgas, T) are related to it. In Table 
[T] we list the adopted parameters for each model. The values at the 
core radius are adopted for the star and gas densities. (We conv ert 
the central densities given in table 1 of Ixsai & Mathew^ 1 19951 to 
the values at the core radius by assuming the radial profiles given 
in that paper.) The stellar mass density p, is used to estimate the 
dust supply rate from AGB stars (Section [2.2. It . and the gas den¬ 
sity Pgas and gas temperature T are used for the initial conditions. 

We start with no dust (Do = 0, where Do is the initial 
dust-to-gas ratio), referred to as the dust-free initial condition, or 
Do = 0.0075, referred to as the dusty initial condition. The latter 
value of Do is typical of the solar-metallicity ISM and the same 
functional form for the grain size distribution as equation (ED [i.e. 
n(a, f = 0) oc fl^®-®] is assumed with the normalization determined 
by equation (0. We expect that the realistic situation lies between 
these two extreme conditions for the initial condition, but that the 
dust-free initial condition is more realistic considering that the dust 
is rapidly destroyed by sputtering as shown below. Thus, we fo¬ 
cus on the dust-free initial condition as a ‘fiducial’ model. As we 
explained above, we investigate two cases for the grain size distri¬ 
bution of AGB dust: MRN and lognormal. 


3 RESULTS 

3.1 Thermal evolution and dust-to-gas ratio 

The thermal evolution of the hot gas is of fundamental importance 
in determining the fate of dust. On the other hand, the temperature 
of hot gas could also be affected by dust cooling. For dust cooling, 
the dust abundance is the key factor. Thus, we first show the time 
evolutions of gas temperature (T) and dust-to-gas ratio (D) in Figs. 
[T]and[3 respectively. 

We compare the results with the different two initial conditions 
for the dust-to-gas ratio. Do = 0 and 0.0075 (dust-free and dusty 
initial conditions, respectively). We adopt the MRN grain size dis¬ 
tribution for AGB dust. As shown in Fig.[D the temperature drops 
more rapidly for the dusty initial condition than for the dust-free 
initial condition because of more efficient dust cooling. The time 
at which the temperature drops down to 10® K, fs, is by a factor 
of ~ 1.5 shorter for the dusty initial condition than for the dust-free 
initial condition. In Table |3 we list ts for each galaxy model with 
different Do and grain size distributions of AGB dust. 
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t [yr] 



t [yr] 


Figure 1. Time evolution of gas temperature T. Note that the gas density 
is always inversely proportional to T during the entire evolution. The upper 
thicker (lower thinner) solid, dotted, and dashed lines present the results 
for Models a, b, and c, respectively, with the dust-free initial condition, 
= 0 (with the dusty initial condition, Do = 0.0075). The MRN grain 
size distribution is adopted for AGB dust. The evolutionary tracks of the 
cases without cooling and with the lognormal size distribution of AGB dust 
for Do both overlap with the thick lines (thus, they are not shown in the 
figure), which means that dust cooling does not affect the thermal evolution 
of the hot gas for the dust-free initial condition. 


Table 2. Time (^ 5 ) at which the temperature drops down to 10^ K. 


IC" 

AGB dust'’ 

Model a 

B (yr) 

Model b 
ts (yr) 

Model c 
ts (yr) 

dust-free 

MRN 

4.9 X 10 '' 

8.3 X 10*^ 

7.3 X 10^ 

dusty 

MRN 

3.3 X 10 '' 

5.3 X lO*^ 

4.5 X 10^ 

dust-free 

lognormal 

4.9 X 10 '' 

8.3 X lO*^ 

7.4 X 10^ 

dust-free 

no dust" 

5.0 X 10 '' 

8.6 X 10*^ 

7.4 X 10^ 


"The initial condition (IC) is chosen from the dust-free and the dusty cases 
with the initial dust-to-gas ratio D = 0 and 0.0075, respectively. For the 
dusty initial condition, the grain size distribution is assumed be the MRN 
(i.e. power law with an index of -3.5; see the text for details). 

^The grain size distribution of AGB dust is chosen from MRN or lognormal. 
"This case is for no dust cooling (i.e. only gas cooling). 

In Fig. we show the evolution of dust-to-gas ratio, D, cal¬ 
culated by equation First we explain the evolution of O for the 
dust-free initial condition. At the beginning of the evolution, the 
dust-to-gas ratio increases because of the dust supplied from the 
AGB stars. The increase of D is saturated because of the destruc¬ 
tion by sputtering. At this stage, the dust-to-gas ratio approaches to 
the value determined by the balance between dust supply and de¬ 
struction as we show later. Indeed, the sputtering time-scale given 
by equation d (note that a < 0.25 pm for the MRN grain size 
distribution) is much shorter than fs, which allows the grain size 
distribution to approach the equilibrium between supply and de¬ 
struction. At the end, the dust-to-gas ratio decreases because the 
compression of the gas due to the temperature drop by cooling en¬ 
hances the sputtering rate (but the temperature is still higher than 
~ 2 X 10^ K). Indeed, the decrease of dust-to-gas ratio coincides 
with the temperature drop (Fig.[TJ. As shown in equation ( fBl l. the 



t [yr] 


Figure 2. Time evolution of dust-to-gas ratio O. (a) The MRN grain size 
distribution is adopted for AGB dust. The solid, dotted, and dashed lines 
present the results for Models a, b, and c, respectively, with the dust-free 
and dust-rich initial conditions (ICs) (lines with enhanced and medium 
thickness, respectively) with cooling included, and with the dust-free IC 
without cooling (thin lines). In the cases with cooling, the calculation is 
stopped when the gas temperature reaches 10^ K (i.e. at time ( 5 ). (b) Com¬ 
parison between the two different grain size distributions of AGB dust: log¬ 
normal (thick lines) and MRN (thin lines). The solid, dotted, and dashed 
lines present the results for Models a, b, and c, respectively. 


sputtering efficiency scales as oc at any gas temperature, while it 
is not sensitive to the gas temperature as long asT > T = 2x 10^ K. 
Because the effect of density is larger than that of temperature, 
cooling event ually enhances the dus t destr uction, contrary to the 
expectation bv iMathews & Brighentil (l2003b . 

In order to show the level of D achieved by the equilibrium 
between dust supply and destruction under the initial gas density, 
we also show the cases where we neglect cooling (i.e. the gas tem¬ 
perature and density are constant) in Fig. [2] In these cases with¬ 
out cooling, the dust-to-gas ratio is the same as the above case at 
the beginning of evolution before the gas cools significantly, and 
it approaches to the value determined by the balance between dust 
supply and sputtering. The dust-to-gas ratio achieved by the equi¬ 
librium is the highest in Model c, primarily because of the highest 
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stellar mass per gas mass (i.e. the highest dust supply rate per gas 
mass). 

We also present the results for the dust-rich initial conditions 
(Do = 0.0075) in Fig. |^. In Models a and b, the hnal dust-to¬ 
gas ratio at is almost the same regardless of Do (note that fs is 
different, though, between the different initial conditions; TableO. 
On the other hand, the effect of initial condition on the final dust- 
to-gas rafio remains in Model c, because the evolution starts with 
lower temperature so the effect of sputtering is smaller than Models 
a and b. Yet, even in the dusty initial condition of Model c, the dust- 
to-gas ratio att = is much lower than the initial dust-to-gas ratio 
(Do = 0.0075). 

Although the large fraction of the initial dust is destroyed in 
any case, the initial condition has an imprint on the temperature 
evolution as shown in Fig.[T] Therefore, if the initial dust-to-gas ra¬ 
tio is as large as the value in the Milky Way ISM, dust cooling has a 
significant influence on the thermal evolution of the gas. However, 
such a high dust-to-gas ratio is never achieved if the dust enrich¬ 
ment and destruction take place simultaneously. 

In order to show the effect of AGB dust, we also examine the 
cases with and without dust enrichment for the dust-free initial con¬ 
dition, finding that dust cooling by AGB dust has little influence 
on the thermal evolution of the hot ISM. Indeed, the evolutionary 
track is indistinguishable between the cases with and without dust 
enrichment (thus, the case without dust enrichment is not shown 
in Fig. [TJ. This is also clear from Table [2] (compare the first and 
last rows): is only slightly longer with no dust cooling than with 
dust cooling. We also confirmed that change of the grain size dis¬ 
tribution of AGB dust to the lognormal size distribution does not 
alter the thermal history. Indeed as shown in Table [T] is almost 
identical between the MRN and lognormal AGB dust. 

In Fig.[2j3, we examine the effects of the grain size distribution 
of AGB dust on the evolution of D. Since the lognormal grain size 
distribution is more biased to large sizes than the MRN grain size 
distribution, it has a longer sputtering time-scale Tsp^ (eauation ll4b . 
The final dusf-to-gas rafio at r fs in the lognormal case is about 
4 times as large as in the MRN case. Nevertheless, since the dust 
abundance is still kept much lower than 0.01 by dust destruction, 
the temperature evolution is not affected by dust cooling even in the 
lognormal case as mentioned above. 


3.2 Grain size distribution 

Now we examine the evolution of grain size distribution in detail. In 
Fig.[3l we show the evolution of grain size distribution for the dust- 
free initial condition. The MRN grain size distribution is adopted 
for AGB dust. In the initial stage at t <K fs, when the change of gas 
density due to cooling is negligible, the grain abundance gradually 
increases, converging to the equilibrium grain size distribution de¬ 
termined by the balance between the supply from AGB stars and the 
destruction by sputtering, and the slope is described by n(a) oc 
(AppendixlBll. We have already seen above that the dust-to-gas ratio 
approaches to the equilibrium value. As is also observed above in 
Fig.[2l the grain abundance decreases at the final stage of evolution 
around t ~ because the density increase by cooling enhances the 
sputtering rate. This is the reason why the grain size distributions 
shift downwards at the final stage in all the three models. 

In Fig. [4] we compare the grain size distributions at t = for 
the three models with the dust-free and dust-rich initial conditions 
(the MRN grain size distribution is adopted for AGB dust). As is 
consistent with Fig. the dust abundance per hydrogen atom (or 
the dust-to-gas ratio) is the largest in Model c: as explained above. 
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Figure 3. Grain size distribution per hydrogen atom for different times. The 
dust-free initial condition is adopted. Grain size distributions are multiplied 
by a'^ to show the grain mass distribution per logarithmic grain radius. Pan¬ 
els (a), (b), and (c) show the results for Models a, b, and c, respectively. 
The correspondence of the line species to the time is shown in each panel. 
The dotted line marked as ‘MRN’ is the MRN grain size distribution with 
a dust-to-gas ratio of 0.0075, shown as a reference for the grain size dis¬ 
tribution appropriate for the Milky Way ISM. Note that the final time (^ 5 ) 
in each model is determined by the time when the gas temperature reaches 
T = 10^ K. The dot-dashed line shows the equilibrium grain size distribu¬ 
tion without dust cooling. 
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Figure 4. Grain size distributions at / = ^5 (after cooling) per hydrogen atom 
(multiplied by a^'). Solid and dashed lines show the results for the dust-free 
(Do = 0) and dust-rich (Do = 0.0075) initial conditions (IC), respectively. 
The dotted line marked as ‘MRN’ is the MRN grain size distribution with a 
dust-to-gas ratio of 0.0075, shown as a reference for the grain size distribu¬ 
tion appropriate for the Milky Way ISM. 


this is due to the highest stellar mass per gas mass (Tabled, that 
is, the highest dust supply rate per gas mass. Moreover, the com¬ 
parison between the cases of different Do in Models a and b shows 
almost no difference, which indicates that almost all the dust that 
existed initially in the hot gas has been destroyed. In Model c, the 
effect of initial condition remains as explained in Section lTTI fFig. 
[He This is because of the low temperature, i.e. low sputtering rate. 
The difference in the initial condition affects at the largest grain 
sizes because the largest grains have the longest sputtering time- 
scale (eguation llTI ). However, even in Model c, the grain abundance 
decreases far below the initial MRN size distribution (shown by the 
dotted line in Fig. II. Therefore, in all cases, the final grain size 
distribution deviates significantly from the initial grain size distri¬ 
bution. 

We also examine the case of the lognormal size distribution 
for AGB dust (equation m. We adopt the dust-free initial con¬ 
dition. The resulting grain size distributions are shown in Fig. 
and their evolutionary behaviour is similar to the case of the MRN 
grain size distribution in the following aspects: The grain abun¬ 
dance increases before significant cooling occurs (i.e. at ? <s: t^), 
converging to the equilibrium grain size distribution (Appendix iBt. 
while it decreases after cooling (i.e. around t ~ t^). Compared with 
the lognormal function, the tail toward small grain sizes becomes 
prominent in the final grain size distributions at t = t^. This tail is 
due to the continuous destruction of large grains by sputtering and 
described by n(a) = constant 1 Appendix iBll. The grain abundance 
is reduced around t = ts because the increase of gas density by 
cooling enhances the sputtering rate, as explained above. 


4 PROCESSING AFTER COOLING 

In the above, we have traced the dust evolution in the hot gas up 
to the point where the gas cools down to 10^ K. After that, sput¬ 
tering does not act as a dust processing mechanism any more, but 
the dust may still be processed by other mechanisms than sputter¬ 
ing. In particular, many authors have argued the importance of dust 
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Figure 5. Same as Fig. but for the lognormal size distribution of AGB 
dust. Panels (a), (b), and (c) show the results for Models a, b, and c, respec¬ 
tively. The correspondence of the line species to the time is shown in each 
panel. Note that the final time (( 5 ) in each model is determined by the time 
when the gas temperature reaches T = 10^ K. The dot-dashed line shows 
the equilibrium grain size distribution without dust cooling. 
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growth by the accretion of gas-p hase metals f or the total dust bud¬ 
get in star-forming galaxies te.g. lDw^ll998h . Since this accretion 
mechanism works in cold and dense clouds, it is worth investigating 
a possibility of dust growth also in the cold gas in elliptical galax¬ 
ies. As mentioned in the Introduction (and shown in Section O, 
it has been argued that the observed dust abundance in elliptical 
galaxies cannot be explained by the internal dust supply from AGB 
stars because of rapid destruction by sputtering. Thus, we exam¬ 
ine dust growth as a p otential additional source of dust in ellip¬ 
tical galaxies (see also [p abian et ah 1994 ^ Donahudl 19951 : 


Martini, Dicken, & Storchi-BergmannI 2013h. 

Under a high pressure in the inner part of elliptical galaxy, a 
cold phase with T ~ 10^ K is prefer red rather than a warm phase 
with r ~ 10'* K dWolfire et al.ll 19951) . Such a cold dense medium 
is indeed suitable for the site of dust growth. Since accretion and 
coagulation (grain-grain sticking) both occur in such a cold dense 
medium, we treat them simultaneously l lHirashitall2012h : however, 
as shown later, coagulation only plays a minor role. Below we cal¬ 
culate the dust evolution driven by accretion and coagulation in the 
cold gas, starting from the grain size distributions achieved after 
cooling of the hot gas. 


4.1 Formulation of accretion and coagulation 


We calculate the evolution of grain size distribution by a ccretion 
and coagulation using the formulation in iHirashlt^ l l2012ll . which 
has a wide application to, for ex ample, polarization predictions 
dVoshchinnikov & Hirashitj|2014h . Again, the silicate dust prop¬ 
erties are assumed, but adopting carbonaceous dust instead does 
not change the evolution of grain size distribution significantly 
dHirashitall2012h . If we consider the pressure equilibrium with the 
hot gas, «H > 10^ cm“^ is achieved for all models at a temper¬ 
ature appropriate for the cold ISM (T ~ 10“ K). Thus, we con¬ 
servatively assume that = 10^ cm“^, but the results for other 
densities can be easily obtained by noting that the time-scale of 
grain processing is simply proportional to . The gas temper¬ 
ature, which determines the thermal speed of the accreting ma¬ 
terial is assume to be 10^ K. The available gas-phase metals (in 
our case, silicon) is esti mated based on the solar abundance of Si, 
(Si/H) = 3.55 X 10-^ dDappeidl200(ih . We also assume that the 
mass fraction of Si in silicate is 0.166. The grain velocity, used to 
estimate the grain-grain collision rate for coagulation, is treated 
as a function of grain radius, and is t aken from the molecular 
cloud case in lYan. Lazarian. & Drainel d2004l) . which has a sim¬ 
ilar gas density to the one considered here. We also impose the 
same c oagulation threshold velocity as adopted in Hirashita_&Jfa3 
(I 2 OO 9 I) (originally based on liChokshi, Tielens. & Hollenbachll 19931 : 
IPominik & Tielensll997h : only if the relative velocity is below this 
threshold, which varies with the sizes of colliding grains, the grains 
coagulate. However, detailed treatment of coagulation does not af¬ 
fect the results, because, as shown later, the abundant gas-phase 
metals makes the role of accretion stronger than that of coagula¬ 
tion in changing the grain size distribution. For the initial condition 
of each model, we use the grain size distribution achieved in the 
hot gas after cooling. More specifically, the initial grain size distri¬ 
bution is produced by multiplying the grain size distributions per 
hydrogen at f = ts, n(a, t^) I tin, in Section[3by uh = 10^ cm“^. The 
time is newly measured from the onset of grain growth in the cold 
gas; to avoid the confusion, we use a different notation for the time, 
t', measured from the onset of grain growth. 



t' [yr] 


Figure 6. Time evolution of dust-to-gas ratio !D by dust growth. The time 
t' is the elapsed time from the onset of dust growth in cooled gas. The thick 
solid, dotted, and dashed lines present the results for Models a, b, and c in 
the case where we use the result of the MRN AGB dust calculation at ^5 for 
the initial condition (i.e. the final state of each model in Fig.f^is used for 
the initial condition). The thin lines with the same line species represent the 
results for the lognormal AGB dust calculation at ts (i.e. the final state of 
each model in Fig.[5]is used for the initial condition). 


4.2 Results 

First, we show in Fig. [6| the evolution of dust-to-gas ratio D esti¬ 
mated by equation at t' < 10* yr (as shown later, the observed 
extinction curves are well covered by the grain size distributions at 
t' < 10* yr). Note that accretion changes O while coagulation does 
not. If we start with the grain size distribution at fs for the MRN 
AGB dust, D begins to increase drastically from f' ~ a few x 10® yr 
and reaches 10“* in a few xlO^ yr. Therefore, if the lifetime of the 
cold gas is longer than 10^ yr, there is a possibility of dust ‘refor¬ 
mation’ using the remnants of sputtered dust grains, which explains 
the observed level of dust abundance as we discuss in the next sec¬ 
tion. The increase of D saturates most quickly in Model c simply 
because the initial dust abundance is the largest, which also means 
that the abundance of available gas-phase metals is the lowest. 

For the calculation adopting the lognormal AGB dust case for 
the initial condition, as shown in Fig.|6] the increase of dust-to-gas 
ratio is not so drastic as in the case of MRN AGB dust models be¬ 
cause the grain size distribution is biased to large sizes (the accre - 
tion time-scale is proportional to the grain radius; lHirashitall2012h . 
This indicates that the initial grain size distribution in the cold gas, 
which is the final grain size distribution of the cooled hot gas, is of 
fundamental importance in determining the efficiency of dust mass 
increase in the cold gas. 

Next, we investigate the evolution of grain size distribution by 
accretion and coagulation. In Fig.|7] we show the evolution of grain 
size distribution, adopting the MRN AGB dust case for the initial 
condition. We observe that accretion continue to increase the dust 
abundance over ~ 10* yr, which is consistent with the evolution of 
O shown above. Because accretion saturates most quickly in Model 
c as explained above, the grain radii achieved by the growth are the 
smallest in Model c. 

The most important feature of grain growth by accretion is that 
the effect appears most significantly at the smallest grain sizes. As 
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Figure 7. Evolution of grain size distribution by grain growth, starting from 
the final grain size distributions in Fig. [3] Panels (a), (b), and (c) show the 
results for Models a, b, and c, respectively. The solid, dotted, dashed, dot- 
dashed, dot-dot-dot-dashed lines present the grain size distributions at f = 
0, 10®, 10^, 3 X 10^ and 10* yr, respectively. We assume «h = 10* cm“*, 
but the time-scale just scales as oc . 


explained in iHirashitJ (l2012h . accretion enhances the dust abun¬ 
dance at the smallest grain sizes, which have the shortest accretion 
time-scale because of the largest grain surface-to-volume ratio. In 
order to isolate the effect of accretion, we also examined the grain 
size distribution with only accretion (i.e. without coagulation). We 
confirmed that there is almost no difference between the cases with 
and without coagulation (we do not plot the results because the 
grain size distributions with and without coagulation are almost in¬ 
distinguishable). This is partly because the effect of accretion is 
enhanced by the large amount of available gas-phase metals, partly 
because coagulation does not proceed beyond a ~ 0.1 p.m at which 
grains have larger velocities than the coagulation threshold (this 
means that the creation of larger grains than ~ 0.1 pm only occurs 
by accretion). 

In Fig. we show the same plot as above but we adopt the 
lognormal AGB dust case (i.e. the final grain size distribution in 
each panel of Fig. |5]l for the initial condition. We observe that the 
change of the grain size distribution occurs more slowly than the 
MRN AGB dust case because the grains are biased to larger sizes 
(recall that the accretion time-scale is proportional to the grain ra¬ 
dius). Therefore, grain growth is very sensitive to the size distri¬ 
bution of the grains incorporated in the cold medium. However, as 
shown later, extinction curves calculated for the lognormal AGB 
dust cases are too flat and are inconsistent with the observed ex¬ 
tinction curves (Section [5.2t . 


5 OBSERVATIONAL PROPERTIES AND DISCUSSION 
5.1 Dust abundance 

For comparison of dust abundance, w e adopt ellip t ical g alaxy 
samples observed by Hersch el from ISmith et aP l l2012h and 
Idi Serego Alighieri et aP ( l2013h . The dust-to-gas ratio is estimated 
for each galaxy based on the dust mass derived from the Herschel 
FIR emission and the gas mass estimated by the sum of H i and 
H 2 masses (the latter is from CO luminosity). The sa mple galax¬ 
ies ad opted are listed in Table For the sample in ISmith et aP 
l l2012h . we select 4 galaxies with (i) a morphological type of el¬ 
liptical galaxies (E); (ii) detection of FIR emission (i.e. a success¬ 
ful estimate of dust mass); and (iii) detection of H i mass. We 
use the upper limits for H 2 mass, but the treatment of H 2 mass 
does not affect the results since these upper limits of H 2 mass 
are significantly smaller than the H i mass. For the sample in 
Idi Serego Alighieri et alj ( l2013h , we exclude galaxies with morpho¬ 
logical type of SO and Sa and only select the objects for which both 
H I and H 2 masses are constrained (i.e. detected or with an upper 
limit obtained). As a consequence, we adopt 4 objects from their 
sample. NGC 4374 is common to both samples. The B-band and 
X-ray luminosities are taken from lO’Sullivan, Forbes, & Ponmarj 

1 I 2 OOII) . 

As seen in Table the observed dust-to-gas ratio is larger 
than lO^'^. However, as shown in Fig. this level of dust-to-gas 
ratio cannot be achieved in the hot ISM, especially in Models a 
and b, which are appropriate for most of the B-band luminosities 
of the sample. The situation does not change even if we change the 
grain size distribution of AGB dust to the lognormal distribution 
(Fig. 1^). This confirms previous conclusions derived by various 
authors that the dust destruction in the hot gas phase prevents the 
dust-t o-gas ratio (or the total dust mass) from reach ing the observed 
level dGoudfrooii & de Jondl995l ; |Patil et al.l2007h . and shows that 
this is also true even if we consider the grain size distribution. We 
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Table 3. Elliptical galaxy sample adopted from ISmith et alj i2012l) an d idi Serego Alighieri et alj l2013l) . 


Name 

Other name 

logifi 

(L©) 

logTx 

(ergs“‘) 

log Mdust 
(M©) 

logMni 

(M©) 

log Mh^ 
(M©) 

log 2) 

Ref.“ 

VCC 763 

NGC 4374, M84 

10.57 

40.83 

5.05* 

8.96 

<7.23 

-3.9F 

1,2,3 

HRS 150 

NGC 4406, M86 

10.66 

42.05 

6.63 

7.95 

<7.4 

-1.43“ 

1, 3 

HRS 186 

NGC 4494 

10.62 

<40.10 

5.08 

8.26 

<7.35 

-3.23“ 

1, 3 

HRS 241 

NGC 4636 

10.51 

41.59 

5.06 

9.0 

<7.02 

-3.94“ 

1,3 

VCC 345 

NGC 4261 

10.70 

41.21 

5.81 

<8.45 

<7.70 

> - 2 . 7 H 

2,3 

VCC 1226 

NGC 4472, M49 

10.90 

41.43 

5.49 

<7.90 

<7.26 

> -2.50“’ 

2,3 

VCC 1619 

NGC 4550 

9.72 

39.78 

5.41 

<7.90 

7.20 

> -2.57“' 

2, 3 


“References: li lSmith et alJ l2012h : 2l ldi Serego Alighieri et all l2013h : 3i lO’Sullivan et al.) bOOlh . 

Adopted from Ref. 1. Ref. 2 gives 5.30 with the same mass absorption coefficient. 

‘ The upper limit of H 2 mass is used, but this does not affect the estimated dust-to-gas ratio since the gas mass is dominated 
by HI gas mass. 

“^We use the upper limits of gas mass to derive the lower limits of dust-to-gas ratio. 


have al so shown that cooling does no t help the dust to survive, con¬ 
trary to lMathews & Brighentil l l2003h ’s expectation: Gas compres¬ 
sion caused by cooling enhances the sputtering efficiency, leading 
to an additional decrease of the dust abundance. In other words, the 
effect of cooling further emphasizes the discrepancy between the 
expected and observed dust abundances in elliptical galaxies. 

Some authors argue an external origin of the gas and dust 
based on the evidence that the kinematic properties of the gas are 
decoupled fro m those of the major stellar component in some ellip - 
tical galaxies ( lForbeslll99ll : ICaon, Macchetto, & Pastorizall2Q0(J). 
Since not all the galaxies fit this external-origin scenario l lForbesI 
Il99lh . it is still worth considering how much dust could be ex¬ 
plained by the ‘internal-origin’ scenario in which the dust is pro¬ 
duced by their own stellar populations and/or inherent dust growth 
mechanisms. 

Regarding the internal origin of dust, we have shown that the 
accretion of gas-phase metals in cold clouds increases the dust- 
to-gas ratio to a level of ~ 10“^ in ~ 10^-10* yr. In this sense, 
cooling is important in the dust abundance , not t hrough survival 
in the hot gas (as in iMathews & Brighentil 1^0031 ’ s original idea) 
but through the formation of cold clouds hosting dust growth. 
The time-scale of the dust abundance increas e by dust growth is 
shorter than the feedback time-scal e (> 10* vr: lKavirai et alj201ll : 
[Pellegrini, Ciotti, & Ostrikerll2012h . on which AGN activity heats 
the cooled gas. Therefore, the ‘recovery’ of dust by accretion is a 
plausible scenario of explaining the observed excessive dust abun¬ 
dance in elliptical galaxies. 

5.2 Extinction curves 

The grain size distributions obtained by our calculat i ons ar e tested 
against the observed extinction curves in IPatil et alj l l2007h . Based 
on the calculated grain size distributions above, we theoretically 
predict the corresponding extinction curves. The extinction at 
wavelength A in units of magnitude (A^) is calculated by weight¬ 
ing the extinction efficiency factor Q^xtCo, A) with the grain size 
distribution n(a) as 

Aa = C f n(a)na^Q,x,(a, A)da, (17) 

Jo 

where C is the normalizing constant, which cancels out when we 
consider the ratio of A^ at two wavelengths later. Although we as¬ 
sumed silicate for the grain p roperty, the grai n growth time-scale of 
carbonaceous dust is similar ( lHirashit3l2012h . The steepness of ex¬ 


tinction curve is often discussed in comparison with the Milky Way 
extinction curve. Therefore, we simply assume a mass fraction of 
0.54 : 0.46 for silicate and carbonaceous dust (graphite), which fits 
the Mi lky Way extinction cur ve under the MRN grain size distri - 
bution jHirashita & Yanl20()^ . We refer to jHirashita & Yanl20()^ 
for the adopted parameters for silicate and graphite. However, since 
we are interested in wavelengths longer than 0.4 pm, the extinction 
curve slope is broadly determined by the grain size distribution with 
a minor effect of the material. 

Now we examine the extinction curves in the MRN AGB dust 
case. In Fig. |3 we show the extinction curves for the grain size 
distributions at 1^ in Fig. [3] (i.e. the initial grain size distributions 
before grain growth) and the variation of extinction curves by grain 
growth (time is measured from the onset of grain growth). We show 
the total to selective extinction = A^IE{B -V)- A^HAb - Ay) 
(the wavelengths at the B and V bands are 0.44 and 0.55 pm, re¬ 
spectively), which is sensitive to the extinction curve slope around 
these two optical bands. The observational extinction curves for a 
sample of elliptical galaxies are taken from IPatil et al.l ( l2Q07h : we 
only chose those galaxies that are classified as elliptical galaxies 
and have extinction data at all bands (B, V, R, and I). 

We observe in Fig. that the initial extinction curves before 
grain growth tend to overproduce since they are too flat [i.e. 
the colour excess E(B - V) is small] compared with the observed 
extinction curves. They are also flatt er than the M ilky Way extinc¬ 
tion curve, which has =; 3.1 (e.e. lDrain3l2003h . because of the 
deficiency of small grains (note that the grain size distribution af¬ 
ter sputtering in Fig. [3 is oc while the Milky Way extinction 
curve is fitted by a grain size distribution oc MRN). After 

grain growth by accretion takes place, the extinction curve becomes 
steeper (or R^ becomes smaller), which is consistent with the pre¬ 
vious conclusi on that grain gro wth by accretion steepens the ex¬ 
tinction curve ( lHirashital2012h . This is because, as explained in 
Section l4^ accretion is much more efficient for smaller grains ow¬ 
ing to their larger surface-to-volume ratio. Since R^ is sensitive to 
the slope at B and V bands, the decrease of R^ is significant if the 
abundanc e of grains with a < A jin ^ 0.07 pm (for A - 0.44 pm) 
increases dBohren & Huffmanll983h . This occurs up to 3x 10’ yr in 
Model a and 10* yr in Models b and c (Fig.|7]l. At 10* yr in Model a, 
almost all grains exceed 0.07 pm, leading to a significant increase 
of R^ (or significant flattening of optical extinction curve), as seen 
in Fig.|9}i. 

It is not likely that we observe the extinction of single-aged 
clouds, but it is probable that we see a mixture of clouds with dif- 
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Figure 8. Evolution of grain size distribution by grain growth, starting from 
the final grain size distributions in Fig. [5] Panels (a), (b), and (c) show the 
results for Models a, b, and c, respectively. The solid, dashed, dot-dashed, 
dot-dot-dot-dashed lines present the grain size distributions at t' = 0, 10^, 
3 X 10^ and 10^ yr, respectively. We assume nn = 10^ cm“^, but the time- 
scale just scales as oc . 





Figure 9. Extinction curves normalized to E{B — V). Panels (a), (b), and 
(c) show the results for Models a, b, and c, respectively. The solid curve 
shows the extinction curves at the final stage of gas cooling in Fig. [3] which 
corresponds to the extinction curve at the onset of grain growth in the cold 
medium. The dotted, dashed, dot-dashed, dot-dot-dot lines present the ex¬ 
tinction curves at 10^, 10^, 3 x 10^, and 10^ yr after the onset of grain 
growth. The squares connec ted by the solid li ne are observational data for 
an elliptical galaxy sample in IPatil et al.N200'^ . The crosses show the Milky 
Way extinction curve as a reference. 
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ferent ages at the same time. Considering this, it is interesting to 
point out that the observed extinction curves are just in the range 
explained by the grain growth. Sputtering in the hot gas tends to 
make the extinction curve too flat; the observed extinction curves 
favour subsequent grain processing that could form grains with 
a < 0.07 pm. Grain growth by accretion enhances the abundance 
of such ‘small’ grains. 

The lognormal grain size distribution of AGB dust produces 
completely flat extinction curve at d < 1 pm, even if we consider 
accretion. Since E{B - V) is slightly negative (i.e. Rx is negative 
and is extremely large), the results cannot be displayed. There¬ 
fore, the lognormal model does not fit the observed extinction curve 
unless there is an extra mechanism of small grain production (see 
Section l53t . 


5.3 Shattering? 

iHirashita & YanI (l2009h also showed that the dust grains may be 
shattered in the war m (~ 10‘* K) diff use medium in the Galactic 
environment (see also lYan et alj|2003) . If some fraction of the gas 
remains to be diffuse after cooling of the hot gas and turbulent ve¬ 
locity is as high as > 1 km s“', grains may be shattered. This would 
enhance the abundance of small grains and would steepen the ex¬ 
tinction curve. If the dust is processed by shattering before being 
included in the cold dense medium, the total surface area of the 
dust is enhanced, resulting in more efficient dust growth by accre¬ 
tion and further steepening of the extinction curve. Therefore, the 
presence of shattering would enhance our conclusion that dust pro¬ 
cessing after cooling is a key to understanding the dust abundance 
and the steepness of extinction curves in elliptical galaxies. How¬ 
ever, in a high pressure environment in the central part of elliptical 
galaxies, such a diffuse warm me dium is not realized as an equi¬ 
librium state dWolfire et al.lll99^ . Therefore, a dynamical model 
is necessary to treat this intermediate temperature range ~ 10“* K, 
which is left for the future work. 


5.4 Implication of the internal origin of dust 

We have shown that accretion in the cold medium is a viable mech¬ 
anism of solving the discrepancy between the expected dust abun¬ 
dance after sputtering and the observed one in elliptical galaxies. 
Moreover, the steepening of extinction curves by accretion in the 
cooled medium is favoured in terms of the steepness of observed 
extinction curves. Although we do not intend to argue that all the 
dust in elliptical galaxies are of internal origin, it is worth noting 
that more careful consideration is necessary before resorting to the 
external origin of elliptical galaxy dust. 

Even if accretion is the most dominant dust source, the dust 
supplied by AGB stars works as seeds for accretion. Therefore, the 
dust supply by AGB stars is still important in obtaining the thor¬ 
ough understanding of the lifecycle of dust in elliptical galaxies. 
The present-day mass loss rate of AGB stars in elliptical galaxies 
can be co nstrained by the m i d-infrared emis s ion em itted by AGB 
dust (e.g. lAthev et ^l2002h . IVillaume et al.l (1201 5l) demonstrated 
that circumstellar dust around AGB stars can account for the mid- 
infrared (8-24 pm) flux in some early-type galaxies. Spectral syn¬ 
thesis models incl uding mid-infrared emi ssion of AGB dust such as 
those presented in ivillaume et al.l ( 1201 5h can be used to constrain 
the dust yield of a population of AGB stars in a galaxy using ob¬ 
served dust emission. This is a crucial component for understand¬ 
ing whether an internal origin is a plausible explanation for the ob¬ 
served dust in early-type galaxies. 


Theoretical dust yields of AGB stars are also tested 
by observations of nearby galaxies. Because of their prox¬ 
imity, the Magellanic Clouds provide us with an opportu¬ 
nity of detailed observation al constraint on the dust forma¬ 
tion i n individual AGB stars ( ^rinivasan_et_^200^Bowrrt^^ 
201 ll; IS riniv asan. Sargent. & l^ ixner 2011 riRiebel et alj 1201 a 



_ Izhukovska & Hennind ( 2013h and 
Schneider et al.l ( l2014h . the current model s of dust production in 
AGB stars used or discussed in this paper Jperrarotti & Gailll^Od ; 
Izhukovska et al.|[2008l ; IVentura et al.ll2014li are broadly consistent 
with the observed dust produc tion rate by AGB stars in the Large 
Magellanic Cloud. However, ISchneider et alj ( l2014h also men¬ 
tioned that there is still a discrepancy between the model and obser¬ 
vation for the Small Magellanic Cloud and that the metallicity de¬ 
pendence of the AGB dust yield is large. Since dust yields of solar- 
metallicity AGB stars have not been tested against observational 
data, observational studies on the dust formation by AGB stars in 
elliptical galaxies, which are, broadly speaking, solar-metallicity 
objects, give an important constraint on theoretical dust production 
models of AGB stars. The understanding of AGB dust yield is also 
crucial even for high-redshift gala xies, in which AGB stars may 
contr i bute to the dust en richment dValiante et al.ll200^ : iGall et al.l 
|201l|;|PiDino etalJzOllh . 


6 CONCLUSION 

We have reconsidered the origin of dust in elliptical galaxies, focus¬ 
ing on the internal origin; that is, the dust originates from the pro¬ 
duction by asymptotic giant branch (AGB) stars within the galaxy. 
The evolution of grain size distribution is consistently solved along 
with the evolution of grain abundance by dust supply from AGB 
stars and dust destruction (sputtering) in the hot interstellar medium 
(ISM). We have confirmed that sputtering is so efficient that the 
dust abundance observed in the far infrared cannot be explained. In 
addition, we have shown that cooling does not help to protect the 
dust from sputtering; rather, gas compression induced by cooling 
raises the sputtering rate so that the dust-to-gas ratio is decreased 
by cooling. Because of the rapid dust destruction, dust cooling has 
little influence on the thermal history of the cooling hot gas, unless 
the hot gas initially has a dust-to-gas ratio as high as ~0.01 (i.e. 
comparable to the value in the Galactic cold ISM). 

In the latter part of this paper, we have considered grain 
growth by the accretion of gas-phase metals in the cold clouds 
formed as a result of cooling of the hot gas, in order to exam¬ 
ine a possibility that the dust abundance is increased or ‘recov¬ 
ered’ by this process. We find that, if accretion lasts longer than 
10 ^/(«h/10 ^ cm“'*) yr, which is shorter than the heating time-scale 
of active galactic nucleus feedback, the dust-to-gas ratio can in¬ 
crease up to ~ 10“^. Therefore, we do not necessarily need to resort 
to the external origin for the observed excessive dust abundance in 
elliptical galaxies, and we suggest dust growth by accretion in the 
cooled gas as a plausible explanation for the dust abundance. Since 
sputtering in the hot gas destroys small grains more efficiently than 
large grains, extinction curves after sputtering are too flat to explain 
the observed curves. The observed extinction curves are better ex¬ 
plained by considering the effect of accretion, which enhances the 
abundance of small (a < 0.07 pm) grains. 

In summary, we conclude (i) that cooling does not protect the 
dust grains against sputtering, (ii) that dust cooling does not change 
the thermal history of the hot gas except for the case in which the 
hot gas initially has a dust-to-gas ratio as large as ~ 0.01, (hi) that 
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the formation of cold dense gas as a result of cooling helps dust 
mass to increase through dust growth by the accretion of gas-phase 
metals, and (iv) that the observed extinction curves are consistent 
with this accretion scenario. 
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APPENDIX A: EVOLUTION OF GRAIN SIZE 
DISTRIBUTION UNDER A VARYING BACKGROUND 
GAS DENSITY 


In the main text, we only concentrate on the single fluid element. 
If we explicitly treat the spatial variation of grain size distribution 
by denoting it as n{a, x, f), where x is the spatial coordinate, the 
continuity equ ation corresponding to equation (TJ is written as (see 
equation 37 in iTsai & Mathewslll995h 


I dn\ 9 

^ 

where u is the velocity of the fluid, and (dldt)a,x is the partial 
derivative for time with fixed a and x. The gas density follows the 
continuity equation written in Lagrangian (comoving) form as 

^+Pga.V-« = 0. (A2) 

Eliminating V • u from equations jAlI l and IA2l l. we obtain 


dn ) 1 dpg, 


Pga; 


dt 


+ -^(an) = S. 
da 


(A3) 


The first two terms are combined to obtain the Lagrangian deriva¬ 
tive (with a fixed): 


1 dp 


9 

-I- —(an) ■■ 
oa 


Pgas d? 

Thus, we obtain equation Q- 


5. 


(A4) 


APPENDIX B: EQUILIBRIUM GRAIN SIZE 
DISTRIBUTION 

As seen in the text, the sputtering time-scale is often much shorter 
than the cooling time-scale. Therefore, it is useful to discuss the 
equilibrium grain size distribution achieved by the balance between 
sputtering and supply from AGB stars (we neglect cooling here). 
Under this equilibrium condition, equation 0 is reduced to 

^[an(a)]^S(a), (Bl) 

oa 

where we set all the time derivative terms as zero. Considering that 
lim„_,ce n(a) = 0, the equilibrium grain size distribution is analyti¬ 
cally written as 

„(a) = ^J‘ S(a)da. (B2) 

If we adopt the MRN grain size distribution feauation ll lb . we ob¬ 
tain 

n(a) = ^^{a-^-^-a-J^), (B3) 

where Omax = 0-25 pm. Here, we have used equation 10 for 5, 
and we can use equation ( 113b for d (note that a is negative). If a is 
significantly smaller than a^ax^ a(a) oc The grain size distri¬ 
bution is less steep than that of AGB dust, since smaller grains are 
more easily destroyed than larger grains. 

In Fig. HD we show the above analytic solution. For compar¬ 
ison, we also show the numerical solution of equation O without 
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Figure Bl. Equilibrium grain size distributions for Models a, b, and c 
(lower, middle and upper lines, respectively). We adopt the MRN grain size 
distribution for AGB dust and neglect cooling (i.e. the gas temperature and 
density are constant). The solid and dashed lines show the analytic and nu¬ 
merical solutions, respectively. For the numerical solution, we plot the grain 
size distributions at t = 10^ yr, which is much longer than the sputtering 
time-scale. 


cooling at t = 10^ yr. Since the time-scales of sputtering and dust 
supply are much shorter than 10^ yr, we expect that the grain size 
distribution has achieved the equilibrium at t = 10^ yr. The ana¬ 
lytic and numerical solutions are almost identical, which confirms 
that our numerical scheme correctly solves the evolution of grain 
size distribution by dust supply and sputtering. Moreover, as ex¬ 
pected above, the slope of the equilibrium grain size distribution is 
described by a power law with an index of -2.5. 

We also show the equilibrium solution for the lognormal grain 
size distribution of AGB dust in Fig. IB2I The numerical solution 
is produced in the same way as above; that is, we used the re¬ 
sult at t = 10^ yr with cooling neglected. The ‘analytical’ solu¬ 
tion is obtained by numerically integrate equation l lB2b . We confirm 
that the numerical results reproduce the analytical results well. At 
a <K flo = 0.1 pm (the central radius of the lognormal distribution), 
n(a) = constant since S (a) drops exponentially in that grain radius 
range. We can confirm this in Fig. lB2l 

This paper has been typeset from a TgX/ FTgX file prepared by the 
author. 
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Figure B2. Same as Fig. lBlI but for the lognormal grain size distribution for 
AGB dust. 
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